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Abstract 

We consider an individual-based spatially structured population for Darwinian 
evolution in an asexual population. The individuals move randomly on a bounded 
continuous space according to a reflected brownian motion. The dynamics involves 
also a birth rate, a density-dependent logistic death rate and a probability of mu¬ 
tation at each birth event. We study the convergence of the microscopic process 
when the population size grows to -1-00 and the mutation probability decreases to 0. 

We prove a convergence towards a jump process that jumps in the infinite dimen¬ 
sional space of the stable spatial distributions. The proof requires specific studies 
of the microscopic model. First, we examine the large deviation principle around 
the deterministic large population limit of the microscopic process. Then, we And 
a lower bound on the exit time of a neighborhood of a stationary spatial distribu¬ 
tion. Finally, we study the extinction time of the branching diffusion processes that 
approximate small size populations. 

Keywords: structured population; birth and death diffusion process; large deviations studies; 

exit time; branching diffusion processes; nonlinear reaction diffusion equations; weak stability; 

Trait Substitution Sequence. 


1 Introduction 

The spatial aspect is an important issue in ecology [301 [H]. The influence of the hetero¬ 
geneity of the environment on the phenotypic evolution has been explored for a long time 
|13[ [T71 |2U] . For example, the emergence of phenotypic clusters under a heterogeneous 
space has been extensively studied dUEZlEI]. In dUEU, the authors suggest that 
clustering and aggregation of individuals can be a consequence of the spatial competi¬ 
tion between individuals. Those phenomena generate a structured population based on 
isolated patches. In m, the authors draw attention to the influence of the boundary of 
the spatial environment. The sensibility to heterogeneously distributed resources is also 
a key point to study the spatial dynamics of population |19| . In this context, the effect 
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of a spatial structure on the evolution of a population is fundamental. 

In this paper, we use a population model that describes the interplay between evo¬ 
lution and spatial structure. We are interested in a macroscopic approximation of the 
microscopic model under three main biological assumptions : rare mutations, large pop¬ 
ulation size and the impossibility of coexistence of two traits for a long time scale. The 
main result of this paper implies a convergence of the microscopic model to a jump pro¬ 
cess that jumps in an infinite dimensional space characterized by the spatial profiles of 
the population. This result is correlated with several works on adaptive dynamics and 
in particular with the model of Trait Substitution Sequence (TSS). Metz and al. |26) 
have introduced this model from an ecological point of view. It describes the evolution of 
phenotypic traits in the case where the ecological time scale, related to the dynamics of 
the population, and the evolutionary time scale, related to the mutations, are separated. 
In the evolution time scale, the model describes the succession of invading phenotypic 
traits as a jump Markov process in the space of phenotypic traits. The link between the 
microscopic model and the TSS model has been completely proved by Champagnat [3] 
in a simpler ecological context. Recently, some papers have generalized this approach 
in the case of an aged-structured population |31| . of a multi-resources chemostat model 
m or of a prey-predator model |6]. But only jSl] deals with some processes with values 
in infinite dimensional spaces, and the age structure is deterministic. In contrast with 
it, we are concerned with the spatial aspect of a population living on a heterogeneous 
environment where the individuals move randomly. 

We study an individual-based model in which any individual birth and death events 
are described. This Markov process has been initially introduced by Champagnat and 
Meleard in [5]. For any time t > 0, each individual i is described by two characteristics 
(XI, Ul). XI represents its location in an open, bounded and convex subset X of with 
a C^-boundary. Ul denotes its phenotypic trait which belongs to a compact subset U of 
The phenotype of an individual does not change during its life time contrary to its 
location. The location space may represent a geographic landscape or a theoretical space 
that describes a gradient of temperature, a gradient of elevation or a resource parameter 
as seed size for a population of birds m In the context of the last example, the phe¬ 
notype may represent the beak size of a bird and it is interesting to study the adequacy 
between the beak size and the seed size when evolution occurs |19) . 

The total population is represented at any time t by the finite measure 
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i=l 

where 5y corresponds to the Dirac measure at y, Nt is the number of individuals at time 
t. The parameter K scales the population size and the biological assumption of large 
population size is stated into mathematics by K tends to -|-oo. 

The dynamics of the process is driven by a birth and death diffusion process, in which 
the motion, birth, mutation and death of each individual depends on its location and 
trait. 

Any individual i with phenotypic trait u moves according to a diffusion process driven 
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by the following stochastic differential eqnation normally reflected at the bonndary dX, 


dXi = V^Id-dBt-n{Xi)dlt ( 1 . 2 ) 

where B is a d-dimensional brownian motion, It is an adapted continnons and non¬ 
decreasing process with Iq = 0, it increases only when XI G dX and the diffnsion coeffi¬ 
cient is a fnnction of the trait. 

We consider a popnlation with asexnal reprodnction. An individnal with location 
X € X and trait u €U gives birth at rate b{x, u). This rate can also be denoted by 6“(x) 
when u is hxed and it is assnmed to be bonnded. The offspring appears at the location 
of its parents. Fnrthermore, a mntation may occnr with probability qkP, making the 
phenotypic trait of the offspring different. The law of the mntant trait is then given by a 
kernel k{x,u, •). The mntation probability p may depend on the trait and the location. 
The parameter qx scales the mntation probability and the biological assnmption of rare 
mntations is stated by qk —>• 0. 

The death rate depends on the characteristics of the individnal and on the com¬ 
petition between all individnals. The natnral death rate is d{x,u). The competition 
exerted by an individnal {y, v) on an individnal (a:, u) depends on the location y and on 
the two traits throngh a competition kernel c : U x X xlA ^ M'*'. For the popnlation 
V = X X]r=i ^ixi,ui) £ Mp{X X U), the competitive pressnre exerted on individnal {x,u) 
is 


1 v-A 1 f 

c-u{x,u) =—'^c{u,Xi,Ui) = — / c{u,y,v)u{dy,dv). 

XxU 

Remark that the competition kernel does depend on y. This spatial dependence yields 
non-trivial mathematical difhcnlties. In Champagnat-Meleard |5], the competition kernel 
depends also on x bnt the long time behavior of the deterministic limit is still nnknown, 
to onr knowledge. Finally, the total death rate is d{x,u) -|- c • v{x,u). As for the birth 
rate, d{x, u) can also be denoted by d“(x) and c(rt, y, v) by c^'^{y) when u and v are hxed. 
Let ns state the assnmptions on the parameters. 

Assumption 1.1. 1. m, h, d, k and c are continuous and non-negative on their do¬ 

mains and b, d and c are Lipschitz functions with respect to x and y. 

2. There exist m, b, b, d, c, c, G M such that for any {x,u,y,v) G (^T x 14)“^, 

0 < < m, b< b{x, u) < b, d{x, u) < d, c< c{u, y, v) < c, k{x, u, v) < k, 

and d is not the zero function. 

3. The sequence of initial measures which belongs to Mp{X xlA), converges 

in law to some deterministic measure denoted by CLnd it satisfies sup^ K[{vq , 1)^] < 
- 1 - 00 . 

4- qx tends to 0 when K tends to -|-oo. 

Before going fnrther, let ns set and recall the notation, which we nse in the entire 
paper. 

Notation 
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• For all X G dX, n{x) denotes the outward normal to the boundary of X at point x. 

• For sufficiently smooth / and for all (x, u) G dX x dnf{x, u) denotes the scalar 
product Vxfix,u) ■ n{x). 

• Cn\x X hi) represents the set of functions / such that / G C^’\X x hi) and 
dnf{x,u) = 0 for all {x,u) G dX x hi. We define Cn^’^{X xhi x [0,T]) similarly. 

• For any / G Cn^’^ {X xhi x [0, T]), fg is the function on X xhi such that fs{x, u) = 
f{x,u,s). 

• For any compact set X, we denote the space of finite measures on X by Mf('K). 

• C^^^(K) denotes the set of all positive Lipschitz-continuous functions / on X 
bounded by 1 and with a Lipschitz constant smaller than 1. 

• We define the Kantorovich-Rubinstein distance on Mf(K) by : for any f, fj. € 
Mf{X), 

u.) = sup / fdv — 

fec^^p{x) Jx 

As X is a compact set, this metric is a metrization of the topology of weak conver¬ 
gence. It is equivalent to the 1^*-Wasserstein distance. 

• B{F,'y) represents the ball of center f and radius 7 in Mf(^) for the previous 
distance. 

• ID)([0, T], Mi 7 ’(X)) denote the space of cadlag functions from [0, T] to Mi?(X), equipped 
with the Skorokhod topology. 

• For any ^ G Mf{X x { m , u }), we identify the two following ways of writing : 

^{dx,dw) = ^'^{dx)6u{dw) + {dx)6v{dw) and ^ G {Mf{X))‘^. 



2 Main theorem 


A full algorithmic description and a mathematical formulation of the model described in 
the previous part are detailed in Champagnat and Meleard [5]. Moreover, a macroscopic 
approximation has been proved as a large population limit. 

Theorem 2.1 (Theorems 4.2 and 4.6 in [5]). Suppose that Assumption \l . 1\ holds. For all 
T > 0, the sequence {f^)k>o of processes belonging to D([0, T], Mi7’(A’ xlA)) converges 
in law to a deterministic and continuous function i.e. ^ G C([0, T], Mp’(A’ x lA)) such 
that supig[o,T]1) < +00 and V/ G Cn^{X x hi), 


{ft,f) = {foJ)+ 



0 XxU 


m^Axf{x, [b{x, u)—d{x, u)—c-fs{x, rt)] f{x, u) ^^s{dx, du)ds. 

( 2 . 1 ) 


Moreover, if hi is finite, for any u € hi and t > 0, f^ti-,u) has a density with respect to 
Lebesgue measure which is a -function. 


The limiting equation (2.1) is a nonlinear nonlocal reaction-diffusion equation defined 
on the space of traits and locations. In mm, the authors have studied the existence of 
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the steady states of similar equations in the context of frequent mutations. Our study 
involves a rare mutations assumption and mutation terms disappear in the limit. The 
stability of the steady states and the long time behavior of the solutions to ( 2 . 1 ) have 
been characterized in | 221 17] in the particular cases of a monomorphic population (all 
individuals have a same phenotype) and a dimorphic population (two traits are involved). 
The stationary states and their stability are described using the following parameters. 

Definition 2.2. For any u gU, we define by, 


= — min 


\h{x) 


/ \V(j){x)\^dx — / — dfi){x)(lP‘{x)dx 

Jx Jx 


( 2 . 2 ) 


where H^{X) is the Sobolev spaee of order 1 on X. is thus the prineipal eigenvalue of 
the operator Ax-+ {b'^-d^)- with Neumann boundary eondition on X. Let G H^{X) 
be the eigenfunction of the previous operator associated with the eigenvalue such that 


[ c^'^iy)riy)dy = HX 
Jx 

According to g^ G C^{X). If g^ > 0, we define the associated measure in Mf{X) 

C'{dx) := g^{x)dx. 


Finally, for any (u, v) G IL, we set 

«-:=^c™(y)r(y)d2/(^r(y)dy) 

As proved in |22) . > 0 is the condition ensuring that a monomorphic population 

with trait u is able to survive for a long time at the ecological time scale. In that case, 
the stationary stable state is described by the positive spatial profile 5 “. 

The dimorphic case implies four distinct stationary states : the trivial state (0,0), two 
monomorphic states and one co-existence state. To ensure that the co-existence state is 
unstable, we set the following assumption. 

Assumption 2.3. Let u be in lA, for almost all v GlA, 


1. either, < 0, 

f > 0 

2. or, < 

1 < 0 . 


Using |22| . we notice that Assumption 2.3 states into mathematics the impossibil¬ 
ity of co-existence of two traits for a long time, this assumption is also known as the 
"Invasion-Implies-Fixation" principle. Under Assumption 2.3, any solution to (2.1) con¬ 
verges either to (^“,0), or to (0,^’'), which are two monomorphic states. More precisely. 
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Condition ensures the stability of the equilibrium (^“,0). Thus, if a mutant popula¬ 
tion with phenotype v is emerging in a monomorphic well-established population with 
phenotype u, it will not be able to survive. Under Condition the authors of 
prove that the deterministic solution to ( 2 . 1 ) converges to the stable equilibrium ( 0 ,^’^ 


whatever the initial condition is. In the light of the previous considerations, we refer 
to as the invasion fitness of the individuals with type u in a resident 

population with type u. Furthermore, the probability of success of such an invasion is 
described precisely by means of the geographical birth position xq of the first individual 
with trait v and the function defined below. That probability is precisely (/>^“(xo). 

Definition 2.4. For any u,v €U, is the function on X such that 

1. If < 0, = 0 for all x£X. 

2. If > 0, is the unique positive solution to the elliptic equation 

irF i^x4>{x) + — (f{x) — J c"^{y)g^{y)dy'^ (f){x) — h"{x)(t){x)‘^ = 0, Vx G X, 

dn(j){x) = 0,Vx G dX. 

(2.3) 

We are now ready to state the main result of this paper. 


Theorem 2.5. We suppose that Assumptions \l.I\ and 2.3 hold. We also assume that the 
scaling parameters satisfy 


KqKlog{K) —> - 1-00 and Kq^e 

K^-\-oo 


KV 


0 , for any U > 0. 


(2.4) 


te[o,r] 


A'^+oo 

converges towards a jump Markov process (Ai)i>[o,r] 


Then for any T > 0^ [^(t/KgK) 

as K ^ -|-oo. At any time t, At belongs to the subspace {f,^6u,u G U} of Mf{X x lA), 
where for any u € U, € Mf{X) is the spatial pattern defined in Definition 2.2 The 

process {At)t>o jumps from the state characterized by the trait u €U to the state charac¬ 
terized by V €U at the infinitesimal rate 

f pb^{x)(IF'^{x)g'^{x)k{x,u.iV)dxdv. 

Jx 

This convergence holds in the sense of convergence of the finite dimensional distributions. 

Remark that describes the spatial distribution of the monomorphic population 
with trait u. The limiting process jumps from a spatial distribution to another one de¬ 
pending on the mutant trait. It models an evolutionary phenomenon using a sequence 
of monomorphic equilibria described by their spatial patterns. 


Although the structure of Theorem 2.5 s proof is similar to the one of Theorem 1 in 
|3], the spatial structure of the process leads us to deal with infinite dimensional pro¬ 
cesses. Two key points of the proof have to be approached differently. The first point 
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concerns the stndy of the process when it is close to a monomorphic deterministic eqni- 
librinm. The aim is to estimate the exit time of a neighborhood of a stationary state to 
(2.11. We give the behavior of the stochastic process aronnd its deterministic eqnilibrinm 
taking into acconnt a small mntant popnlation and the possibility of other imitations. 
Thanks to it, we avoid the comparisons nsed in [3], where the behavior of the resident 
popnlation process is compared with the behavior of a theoretical monomorphic popnla¬ 
tion evolving alone. Those comparisons are mnch more involved when the popnlation is 
spatially strnctnred. Moreover, to estimate this exit time, we have to stndy a large devi¬ 
ation principle of the stochastic process (r'/^)t>o aronnd its deterministic limit (2.1 ) when 
K is large. The large deviations studies for processes combining diffnsion process and 
jnmps are still nnresolved, to onr knowledge. Those stndies have thns their own interest. 
The second point which is approached differently concerns the stndy of small popnlation 
size processes. The aim is to nnderstand the dynamics of a popnlation descended from 
a mntant which has appeared in a well-established monomorphic popnlation. As long 
as the mntant popnlation size is small, the competitive terms between mntants can be 
neglected. Thns, the dynamics of the mntant popnlation can be compared with the dy¬ 
namics of a branching diffnsion process. We describe finely the snrvival probability of a 
branching diffnsion process by means of the eigenparameters defined previonsly and we 
link it with the conditions presented in Assnmption |2.3| 

In Section we explicit the npper bonnd of the large deviation principle by nsing 
ideas in |9l |25l [3T] . Then we stndy the fnnctional rate associated with the large deviation 
principle. Section deals with the exit time of a neighborhood of a stationary state to 
(2.1). In Section]^ we stndy a branching diffnsion process. First, we evalnate its prob¬ 
ability of snrvival. Then, we characterize the scale time nnder which its size is of order 
K. Section is devoted to the proof of Theorem 2.5 We detail two key propositions. 
The first one deals with the dynamics of the individnal-based process in the case where 
only one trait are involved. The second one gives the dynamics of the process after the 
time of the first mntation bnt as long as only at most two traits are involved. Finally, 
Sectionpresents a nnmerical example that illnstrates Theorem 2.5 


3 Exponential deviations results 


In this section, we are concerned by the large deviations from the large popnlation limit 
(2.1) for the process (r'/^)ig[o,T] when K tends to -(-oo and qx tends to 0. First of all. 
Theorem 3.1 gives the npper bonnd of the large deviations principle. This theorem 
involves a rate fnnction. Let ns first explicit it. That requires specific notation which 
will be only nsed in this snbsection : let ns fix T > 0, 


• S = X xUx {1,2}. 

• Ip is the mapping snch that for any fnnction / G xlA x [0,T]), for any 
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(x, u,7r,t) € £ X [0, T] 


f{x,U,t) if TT = 1, 
— f{x, U, t) if TT = 2. 


Tp{f){x,U,t,TT) = 

• For all v = ii^t)te[o,T] ^ ®([0) Mp{X x U)), we define the positive finite measure 

{dx, du, dvr) = [b{x, u)6i{d7T) + {d{x, u) + c - ut-ix, u))S2{d7r)~\ ut-idx, du). 

(3.1) 

• Finally, we introduce the log-Laplace transform p of a centered Poisson distribution 

with parameter 1, p{x) = — x — 1, and its Legendre transform p*, 

p*{y) = {(y + 1) iog(y +1) - y)i{y>-i} + i{y=-i} + oo • 

We are now ready to define the rate function in which we shall be interested : for all 
^0 G Mf{XxU) and € D([0, T], x L/)), 

( sup if t'o = Co 

:= } fGC^’°’\xxUxlo,T]) ( 3 , 2 ) 

I, + oo, otherwise. 


where 


:={ixp, fp) - {pq, / o ) - j {niA^fs + m\V^fs\‘^ + ^,Us)ds 

rT 




When there is no ambiguity, we will write l'^{v) instead of 

Theorem 3.1. Suppose that Assumptions \n\ holds. For all a > 0, G Mf[X x lA), 
for all compact set C C for all measurable subset A o/D([0, T], Mi?(T x U)) 

such that there exists M > 0 with A C {i^| sup^gjQ (r't, 1) < M}, 


limsup ^ sup logP^if(i^^ G A) < — inf (3.3) 

K^+oo di ° ^&C,u&A 


where M|f = {X with IV e N, {xi,Ui) € X xU}. 

Proof. We will show the following upper bound 


limsup — \og¥{u^ G A) < — inf liAv), 
K ..^ /I ‘SLi 


K^+oo 


uGA 


(3.4) 


indeed (3.3) can be directly deduced from this bound by a similar reasoning as in the proof 
of Corollary 5.6.15 in Dembo and Zeitouni [9]. To prove (3.4), we need the exponential 
tightness of the process (r'/^)te[o,T] which is described by the following lemma. 
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Lemma 3.2. Suppose that Assumption \1.1\ holds, and that there exists Cinit > 0 sueh 
that supj^gf!}(z^(f, 1) < Cinit O'-S-- Then for all L > 0, there exists a eompaet subset Cl of 
the Skorohod spaee x Li)) sueh that 

limsup^ logP(i^'^ 0 Cl) < —L 

K^+oo hC 


We do not detail the proof of Lemma 3.2 as it may be easily adapted from laiiHiEi]. 
Then, set = inf{t > 0, 1) > M}. Note that Lemma 3.2 is also true for K)t>o- 


Using a proof similar to Theorem 4.4.2 of we deduce the inequality 


limsup^ logP(i^^ ^ A) = limsup^ logP(i/^ k ^ A) 

K^-\-oo ^ K^-\-oo ^ ^ 

< - inf ( sup - H{lf’'^))'), 

f&Cn’°’^{XxUx[0,T]) ^ 


where H{lf''^) = limsup;^^_,_oo 4 logE[exp(iir/-^’^(z/^ ^)]. Let us show that H{lf''^) = 
0. Let 


A/'t = exp (k) 


- QkK 


r'TArf 


vf,pb^{x){ I (f{fs){x,w)k{x,u,w)dw - (f{fs){x,u] 

u 


ds 


where (j){x) = x + p{x) = e* — 1. Ito’s Formula implies that {Nt^T > 0) is a local 
martingale. The dehnition of implies that A/'t is bounded. So it is a martingale of 
mean 1 and there exists a constant C(||/||oojAf) > 0 such that 


exp(-(?,^A:C(||/||oo,M)) <E 


< exp(gi^A:C'(||/||oo,M)). 


[exp 

We conclude easily, since qx tends to 0 when AT —>■ +oo. □ 

The aim is now to write the rate function under a non-variational integral formulation 


which is more workable than that of (3.2). Firstly, this integral formulation is convenient 


to use Chasles’ Theorem. Secondly, it will be used to bound from above the distance 


between a solution to (2.1) and any i/, this upper bound is proved below in Proposition 


3.4 Those two points are required to prove the results about the exit time in Section]^ 
Before writing the non-variational formulation, let us dehne two functional spaces. 

• The Orlicz space associated with p* is Lp*^T the set of all bounded and measurable 
functions h on £ x [0, T] such that 


qT := inf > 0, 


p* [-—^ ) dp’^ds < 1 > < -1-00. (3-5) 

fxfo.T] V “ 


The Orlicz space associated with p is dehned on the same way. 
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is the set of functions h € L‘^{X xU x [0,r],M^) such that 




fT \ 1/2 

2{vs-,'m\hs\^)ds\ < oo. 


'0 


(3.6) 


Theorem 3.3. Suppose that Assumption holds. Let T > 0 and v G ®([0, T], Mp{X x 
hi)), such that Ivq{^) < +oo, then there exist two measurable functions {h^, / 12 ) G Lp*,T x 
such that for all f G Cn^'^{X xU x [0, T]), 


'XtJt) = (t'Oj/o) + 

rT 


+ 


10 


/ (1 + hi{x, u, s, 7r))f;{f){x, u, s, 7r)diJ,fds 
JO Js 

[ fm“Aa;/s(x,u) + 2m“/i2(x,u,s) • Vxfs{x,u) + ^^{x,u)]vs{dx,du), 
JxxU ^ OS J 

(3.7) 


and the rate function can be written as follows 


= 


10 


T r rT 

I p*(hf)dfXgds + / m(i/s, |/i 2 p)ds <+ 00 . 
s Jo 


(3.8) 


The proof of Theorem 3.3 uses convex analysis arguments which can be adapted 
from Leonard [21 El]. We do not detail its proof but we give the main ideas. For all 
T G D([0, T], Mi 7 ’(/t’ X hi)), I^^iu) is equal to the Legendre transform F* of 

'(V^,V,)(C2A1(A'x/7 X [0,r])) ^ M 

(91,92) H- fg p(gi)dp’'^ds + f^{ns,m\g2\^)ds^ 

(3.9) 

at a well chosen point l,^. If 1,^ belongs to the interior of the set domT* of linear maps 
I with r*(/) < + 00 , we can exhibit 1^ by means of the derivative of F, the Legendre 
biconjugate of F. Studying directly F is difficult. The key point is thus to work on 
the product space Lp^p x Cff^. In this way, we can study the Legendre biconjugate of 
an extension of F on that space, in order to deal with the diffusive part and the jumps 
part separately. The diffusive part is treated using ideas of Dawson and Gartner | 8 | and 
Fontbona mi whereas the jumps part is treated using ideas of Leonard |23[ I24| . The 
next step is to deduce the Legendre biconjugate of F by restricting the dehnition domain. 
Finally, to deal with points v for which G does not belong to the interior of domV*, we 
use a continuity argument similar to that of Theorem 7.1’s proof in |25| . 


The last result of this part gives an upper bound on the distance between a solution 


to (2.1) and any u, this bound is used in Subsection 4.2 


Proposition 3.4. Let T > 0 and M > 0. There exists C(T,M) such that, for any n 
satisfying sup^^p{ut, 1) < M and for all {f,t)t>o solution to (2.1) with the initial condition 
f,o = ^0, 

sup Wi{ut,it) < C{T,M) (l^{v) + . 

te[o,T] V / 
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Proof. Let u be such that (t't,!) < M for all t G [0,T]. If = 0, i.e. ut = or if 

l'^{v) = + 00 , the result is obvious, so let us assume that 0 < < +oo. Theorem 3.3 

implies the existence of (hi, / 12 ) G Lp*^T x such that 


{^) = [ [ p*ihi)d^gds+ [ m{us,\h2\‘^)ds. 

Jo Je Jo 

We easily deduce that for any t <T, 

||h 2 ||£ 2^1 < II^ 2 ||£ 2 ^' 7 ’ < 2 / (l/). 

Let us also find an upper bound on ||/ii||p*^'r. Note that for all x G M, 

J if a > 1, p*{\x\/a) < p*(|x|)/q; < p*{x)/a, 

[ if 0 < a < 1, p*{\x\/a) < p*(|x|)/q;^ < p*{x)lo?. 


(3.10) 


(3.11) 


(3.12) 


Moreover, the non-variational formulation (3.10) implies that/^(z 2 ) > p*(hi)dpgds. 

Thus, using ( 3.12[ ) and the definition of the norm ||.||p*,T hr (3.5), we obtain that if 
/^(z/) > 1, /g f^p* dp^ds < 1, i.e. ||hi||p*,r < h^(i^), and if < 1, 

fo fsP* (|hi|/A//^(n)^ dpfds < 1, i.e. ||hi||p*,T < Thus, for any t<T, 


||hi||p.,t < ||hi||p*,r < • 


(3.13) 


Let {f,t)t>o he the solution to (2.1) with initial condition vq. We want now evaluate 


Let us denote the semigroup of the reflected diffusion process which is the 


solution to (1.2) with the initial condition x and the diffusion coefficient by {P^)t>o- 
Using Theorem |3.3[ we find the following mild formulation for {i't)t>o in a similar way 
to Lemma 4.5 in m : for all / G C^^P{X x U), 

[ [ 'if{Pt-sf)dPsJs+ f [ V'(-Pt-s/)hi(i/r^ds + f {us,2mV^Pt_J-h2)ds. 

Jo Je Jo Je Jo 

(3.14) 

In addition with a mild equation for (^t )i> 0 ) we deduce that for all / G C^®^(T’ x hi) and 


for all t < T, 


\Wt-Ct,f)\ = 



fjiPt-sDidPsds - dpids) + [ [ 'if{P^_J)hidp'',ds 

Jo Je 


/o Je 

+ [ {us,2mVxPt-sf ■ h2)ds 

Jo 

<Ci f sup >Vl(l/r-,^r)hs + ||V’(-F)-_./)||p,t||hi||p*,t + ||Va;Pt_,/||£2_J|/i2||£2_i. 

Jo rg[0,s] 

(3.15) 

The second line is a consequence of Holder’s inequality (see for example Theorem 6 of 
Chapter 1 in |2^ about Holder’s inequalities in Orlicz spaces). 

Furthermore, the following Lemma insures that P^f G C^'‘P{X). 
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Lemma 3.5 (Part 2 of |32]). As Af is a convex set in for all f G u £ U, 

and t G M+, P/*/ G C^^P{X). 

Let us now find an upper bound on \\'f’{Pi-.f)\\p,t- f belongs to C^^P{X x U) and 
supjg[o^ 7 ’] (t't, 1) < M. So, for all a > 0, t < T 



mPt-sf)\ 


a 


dfigds 


< 



< TM[b + d + cM]p 



so, for any t <T, 


mPt-.f)\\p,t< 


f —— 

\TMlb + d + cM] 



:=C 2 . 


Furthermore, Lemma 3.5 implies also that for any t <T, 


W'^xPi- f\\c 21 = / ‘^{^s,m\VxPt-sf\‘^)ds < 2mMt < 2mMT := C 3 . 
’ Jo 


Using the last two inequalities with (|3.1ip, (|3.13|) and (|3.15|), we find 

ip Wl{Ur,^r) <Ci [ 

re[l 


sup Wl{Ur,^r) <Ci f sup Wl{Vr, ir)ds + C 2 ({u) + ■ 

'e[ 0 ,r] Jo re[0,s] V v y v 


We use Gronwall’s Lemma to conclude. 


□ 


4 Lower bound on the exit time of a neighborhood of the 
stationary state 


In this section, we assume that initially, two traits u and v are involved. The stochastic 
process starts in a state Vq = such that is close to and there exist 

only a few individuals with trait v. Since the considered initial state is close to the 
equilibrium (^“, 0) and according to Theorem 


2.1 


the dynamics of the stochastic process 
is close to the equilibrium (^“,0) on a finite interval time when K is large. Our aim 
is to control the exit time of the stochastic process from a neighborhood of the 

stationary solution in Mp[X) when K is large and is small. We define the exit 

time by : 


for all 7 > 0, P; 


K 


= inf{t>0,Wi(z.f’“,n>7)}- 


(4.1) 


Theorem 


4.1 


gives a lower bound on . The lower bound involves the first time when 


a new mutation occurs and the first time when the n-population size is larger than a 
threshold: 


Sf' = inf{t > 0, 3w 0 {n, n}, I'fiX x {tc}) / 0}, 
for all e > 0, = inf{t > 0, , 1) > e}. 
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(4.2) 

(4.3) 


















Theorem 4.1. Suppose that Assumption 1.1 holds and that > 0. Let 7 > 0 such 
that 7 < , and if H'^k™ — > 0 , 7 satisfies also the assumption 7 < 

I I Then, there exist 7 ' > 0, e > 0, and y > 0 such that, if 


^UU t~VV _ ^vu ^uv 


,K 


= t'n 


K.u 


+ with < 7 ' and < e, then 


lim (Rfi > e^^ ATT ASfi) = 1. 

Thus, a well-established monomorphic population u is minimally affected during the 
emerging of a mutant population v. 

The assumptions on the radius 7 of the neighborhood ensure that there exists only 
one steady state in the neighborhood. 

The result is proved using ideas similar to the ones of Freidlin and Wentzell |15| . 
in our framework, the difficulties come from the continuous space motion. Firstly, our 
processes have values in an infinite dimensional space, thus, the required deterministic 
results are much more involved, see Subsection 4.1 Secondly, we deal with two kind 


of randomness : jump process and spatial diffusion process. The end of the section is 


devoted to the proof of Theorem 4.1 


4.1 Stability for the weak topology 

This subsection deals with the deterministic solution to (|2.1[). We denote by {f,t)t>o 


the solution to equation (2.1) with initial condition ^0 £ Mf{X x {n,n}). In this case, 
ft ^ A1 f{X X {u, n}) for all t > 0. We prove that, as long as the size of the n-population 
density is small, the n-population density stays in a Wi-neighborhood of its equilibrium 




Proposition 4.2. Suppose that Assumption \n\ holds. Let 7 > 0. There exist 7 ' > 0 
and e' > 0 such that for any ^0 = ffidu + 

for all t<tF = mi{t > 0 , 1 ) > e'}, < 7 / 2 - 


The proof of Proposition 4.2 implies two main difficulties. First, using ideas similar 
to Part 3.3 of [22], we can prove that the solution ft to (2.1) stays close to if the initial 


condition admits a density which is close to the density of for the L^-distance. 
However, this is not sufficient since we will deal with discrete measures later. Thus, we 
need to enlarge the result for Wi-distance. Secondly, we are concerned with the trajec¬ 
tories of the M-population process. Even though the n-population size is small, it does 
have an impact on the death rate of individuals u which we cannot ignore. 

The proof is divided into three steps. Firstly, we study how fast a solution with initial 
condition close to moves away from in Wi-distance during a small time interval 
[0,to]- Then, as Iq > 0, ff^ admits a density and so, we can compare the Wi-distance 
and the L^-distance of the densities between fjf and f^. We finally prove a L^-stability 
result similar to the one of Part 3.3 in |22| but including the n-population process with 
a small size. 
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Proof of Proposition 4-2 First we may assume that e' < 1 and 7 ' < 7 . Hence, there 
exists M > 0 such that any considered initial state satisfies (^Qj 1) < 

We fix to > 0 and we start with the first step. On the one hand, we can find an upper 
bound to sup^g[o,t] ) !)• Indeed 


+ b / {Ca)ds, 

Jo 


(4.4) 


and using Gronwall’s Lemma, we deduce that sup^g[o,t](0) H) — for all t > 0. On 

the other hand, using (2.1) with U = {u,v} and a mild formulation similar to Lemma 
4.5 in m, we find that satisfies: for any / G C^^P{X), 

- r, /) =(Co - ^“/) + f\cs - r, -d^- • e)pr-sf)ds 


+ / c^^-{e-cm,ptsf)ds- / {c^'’-c){c,pr-sf)ds. 

Jo Jo 


For any g Lipschitz-continuous, we denote by H^HLip the smallest constant such that 
g/Wghip^ Since sup^gfo^^^,] 1 ) < e', using Lemma 

distance Wi, we obtain that, for all t < tf^i, 


3.5 


and the definition of 


i(er-r,/)i<wi(^o“,n + (5+rf+c(r,i)) / wi{^s,ods 

Jo 

+ \\c\\Lip sup ( [ Wi{C,C)ds + e''] . 

rS[0,t] \Jo / 


(4.5) 


Finally, (4.4), (4.5) and Gronwall’s Lemma imply that there exist Ci,C 2 independent of 
e' and 7 ' such that 

sup Wi(C, C) < (Wi(eo, n + < ( 7 ' + e'C 2 )e^i*°. (4.6) 

re[0,toAt^,] 


According to (4.6), we have to choose 7 ' and e' such that ( 7 ' + e'C 2 )e^^^° < 7 / 2 . Note 
that if for all € Mf{X), t^' < toi the proof of Proposition 4.2 is complete. In what 
follows, let us assume that > to for the considered initial state ^q. 

The next step is to compare the L^-distance and the Wi-distance between and 
According to Theorem 2.1, for any to > 0, has a Lipschitz-continuous density with 
respect to Lebesgue measure on X that we denote by gf^(x). In addition with the fact 
that 5 “ G C'^(A’), we have 


hi - rwi. = [ hiix) - r{x)fdx < mai^emiwup + wrwup). ( 4 . 7 ) 

Jx 

Let us bound \\gf;^ Wup from above. For any f > 0, we define hf{x) = gfix) exp(f*° (c““ • gf + c“'" 
The exponent of the exponential term is positive and independent of x, thus ||i,ip < 


■a)ds). 
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Il^^o \\Lip- Furthermore, according to Part 4 of chapter 5 in [16], (x) = Ftp (x, y)Co{dy) 

where P is the fundamental solution to the system 

' dtV = m“Ar + (6“(x) - d“(x))r on d:” x M+, 

< dn^ = 0 on dX X M“'“, 

_ r(0, dx) = (o{dx). 


As to > 0) l|Fto|| l°°(a’) s-iid IIVPtQ are bounded from above and there exists C 3 

such that 


UoWup < \K\\up < < CsM, 


(4.8) 


where M has been dehned in the beginning of the proof. Combining (4.6), (4.7), 
and the fact that Wi(.^o,f“) < 7 ', we hnd Ci{e','y') > 0 such that 


(4.8) 


\\9l-r\\L^ < ((Wi(eo“,r) + ^^2e')e''^*°(C3M + ||rHi, 

< ((7' + C 2 e')e^i*°(C 3 M + IIl,^)):= ^4(6', 7')- 


We now deal with the last step of the proof. We dehne (A^, Ak)k>i the spectral basis 
for the operator m“Aa,. + ( 6 “ — dA). with Neumann boundary condition such that (Afc)fc>i 
is a non-decreasing sequence with := Ai > A 2 > A 3 ... and {Ak)k>i is an orthonormal 
basis of L^(A’). Note that is equal to ||g'“ 11 ^ 2 ^ 1 . Let us express g'^ — 5 “ in the basis 
{Ak)k>i 

+ 00 

9t{x) = 5 “ (a;) + ^ai{t)Ai{x), 
i=l 


and denote for all t G M"*", 


/3(t) 1 -|- 


ai{t) 

lid" 


IL 2 


From (2.1) and the representation of — 5 “ and dtg^ with respect to the basis {Ak)k>i-, 
we hnd the following dynamical system 


dtakit) = ak{t) I^Afc - c^^gt - - 5 “)^ , for all k>2, 

dtm = m - E ■ 


(4.9) 

(4.10) 


The last step of the proof consists in proving that if e' and 7 ' are sufficiently small such 
that 

maxjl, I • (e'c-h3||c““||i2C'4(e',7')) < 2p^} ’ 

then for all t < Wi( 5 “, 5 “) < 7 / 2 . Let us hx iP > 0 such that 


e'c + 3||c““||i2C4(e',7') <iF/2< (F“-A2)/2, (4.12) 
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and 

tmax = inf{i > io, I [ - 5“)l > K}. 

Jx 

Let us prove that tmax > ie'- At t = to, I/;,:■ - 5 “)! < ||c““||i 2 C' 4 (e', 7 ') <-^72 < AT. 

A continuity argument implies that tmax > ^o- Then, using ( |4.9| ), we deduce that for any 
k>2, for any to < t < tmax A t^/, 

oik{tf < a(to)^ exp(2(Afc - + Ar)(t - to)) < a(to)^. 


Thus 

+00 


A 

^a.(t) / (c”“A) 

i=2 


< 


< 



1/2 


1/2 


(4.13) 

Inserting ( |4.13 ) in Equation (4.10) implies that, for all t < tmax A t^/, 

/3(t) - c'c- C4(e',7')l|c““||L2 - H^m) < 5t/3(t) 

< /3(t) + e'c + C4(e',7')l|c"“||L2 - F"/3(t)) . 

Moreover, |/3(to) —1| || 5 “|| 2,2 < < 74 ( 6 ', 70- Using properties of logistic equations, we deduce 
that for all to < t < tmax A t^/, 

|ai(()l = l/3(*) - lllirili^ < l^(e's + C4(e',7')l|c"“lli.) + C4(e'.7')- (4.14) 


Thus, for all to < t < tmax A t^/, (4.13) and (4.14) imply that 

+O0 

c^^{gt-r) 


IX 


<i«i(t)i^!^ + 




^a,(t) [ (c““A.) 

i=2 




< 6'c-+2C4(6',7')l|c““k^ + ^i^^4(6',7')- 


By dehnition, < ||c““||£ 2 ||g '“||£,2 and (4.12) gives for all t < tmax A V, 


c““(5r-r) 


/A" 


<e'c + 3C'4(e',7')l|c““||L2<^. 


(4.15) 


Inequality (4.15) implies that tmaxf\t^i < tmax, that is, tmaxf^t^i = tf:> and all inequalities 
proved above are true for all t < t^/. In addition, (4.11), (4.13) and (4.14) imply that, 
for all t < t^/, 

wiisr.r) < liar - 9“iii.iiiiii= < iiiiin nai(t)i +1 

7 

2 ' 

That ends the proof. □ 


< ^e'c + 3C,{e',j')\\cn\L^) < A. 
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4.2 Exit time 


This subsection is devoted to the proof of Theorem |4.1| We split the proof into three 
lemmas similar to the ones in Dembo and Zeitouni [^. 


Let 7 > 0 sati sfying the assumptions of Theorem 4.1 We consider e' and 7 ' as in 


Proposition 4.2 and set e = §■ and p = y < 7 . , Si and have been defined by 


■^K 


(4.1), (4.2) and (4.3) and let us define 


T = inf{t > 0, Til- 

Lemma 4.3. Under Assumption we have 


lim lim sup — log 

t-!- + 00 K 


sup .. J,. 

i^^e(S(r,7)xS(0,e))nM|f ° 


k{t A A Si > t) = —00. 


K 


That is, the probability that the process stays a long time in the ring 5(^“, 7 ) \ 
is exponentially small. The proof requires a comparison with the deterministic 


paths of Equation (2.1), the difficulty is to prove that there exists a finite time after 


which all deterministic paths starting in the ring are out of the ring. The fact that the 


probability is exponentially small is a consequence of Theorem 3.1 


Lemma 4.4. Under Assumption [U there exists V > 0 sueh that 


lim sup — log 
K^+00 iA 


sup 


> 7,r < A ) < -E 


,K 




Once again, the proof is based on a comparison with the deterministic paths, it is a 


consequence of Proposition |4.2| and Theorem 3.1 


Lemma 4.5. Under Assumption [m for all C > 0, there exists T{C,p) > 0 sueh that, 


lim sup — log 


K^-\-oo 


K 


sup 


i/|f’“eB(C“,7)nM|f \te[o,T(C,p)] 


sup 


Wi{uf’\u^’^)>p 


< -C. 


This lemma means that there exists a small time interval during which any process stays 
close from its starting point with an exponentially high probability. The stochastic pro¬ 
cess includes jump parts and diffusive parts. Thus, we study not only the size of the 
population process during a small time interval, but also a sum of reflected diffusion 
processes. 


Theorem 4.1 is proved using the two last lemmas 4.4 and 4.5 We do not detail its 


proof as it can be adapted from [9]. The main change is that the proof has to be done 
on {t < A S'f} to ensure that the u-process size , 1) is small and that no other 


mutation appears, but Lemma 4.4 is sufficient to circumvent this difficulty. 
There remains to prove the three lemmas. 
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Proof of Lemma 4-3 Remark that if < p for all RT > 1, the result is obvi¬ 

ous. Otherwise, let us dehne the following set: 

A{T) = adh{v E B{[0,T], MF{Xx{u,v})),yt E [0,T],Wi{uf,C) ^]p,l[ and {uf,l) < 4, 

(4.16) 

where adhS is the closure of the set S. Remark that E A{T)} = {t AT^ A > T} 
a.s. and that the set C = adh{B{fA^l) ^ -®(0,e)) is a compact set of Mf{X x {u,v}) 
as X is bounded. By applying Theorem 3.1 with the closed set A{T) and the initial 
compact set C, we hud 


lim sup — log sup 
K^+oo K 




v(iA{T) 


Thus, the lemma will be proved if we show that infj^g^^)/'^(z/)—>--|-oo as T tends to -|-oo. 
To this aim, we will first show that any solution to ( |2.1[ ) cannot belong to ^(Tq) if Tq is 
large enough. Precisely, we set 6 e] 0, p/2[, and we will prove that there exists Tq > 0, such 


that any {ft)t>o solution to (2.1) with an initial condition satisfying Wi(^q,^“) E]p, 7 ) 
and (^Q, 1 ) < e satishes 

inf sup > 5. (4.17) 

^^XTo)te[0,To] 


Assume that (4.17) holds. Since supj^g_ 4 ( 7 |j) supjg[o,ro] (fo i® bounded, we can use 
Proposition |3.4| to deduce that there exists C > 0 such that for any f E .4,(To), 


we 


6 < C{I^° (f) + y^/T>(z/)). As X 1 — 7 ' X -|- ^/x is a bijective function from to 
hnd a constant C{S) >0 which is a lower bound on I'^°{f) and 7^°(z^) > C{5). 

Finally, for T > Tq and f E A{T)^ we decompose 1 / as a sum of n := \T/Tq\ terms : 
using Chasles’ Theorem and the non-variational formulation (3.8) of > 0, we hnd a 
sequence (i/*)j=i,E .A(To) such that 


i=l 


n —;■ - 1 - 00 . 

T^+oo 


There remains to prove that (4. ^ holds. Let >Vi(^o,^“) E]p, 7 [ and {(,q,1) < e, 
and let {^t)t>o be the solution to (2.1) with initial condition ^0 = ^q^u + ^q^v Using 
Theorems 1.2 and 1.4 in |22| . we know that E Mf(X x {u, x}) converges to a stationary 
state. Theorem 1.2 in |22] implies that this stationary state is either or or a 
state with coexistence but the assumptions on 7 insure that those stationary states do 
not belong to (R(^“, 7 ) \ B{^^,p)) x B{0,e). Therefore, there exists such that 


44^1 7 ^ [7 + 25, p — 25] or , 1) > e -|- 25. 


(4.18) 


However depends on ^q. Thus, we will use a compactness argument to hnd a uniform 
time and conclude. First, using Gronwall’s Lemma, we obtain : for M > (^“, 1)-|- 27 -|- 2 e. 


sup ( 6 , 1 ) < ( 6 ,l)e^^^° < 

4e[o,7e(j] 


(4.19) 
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Then we show that two solutions {^t)t>o s-nd {Ct)t>o to (2.1) which are close initially stay 
close durine a short time. Indeed, usine the mild equation of Lemma 4.5 in [5], for any 
/ e C^^P{X X {rt,?;}) 

i^t-CtJ) = {^o-Co,Ptf)+ [\^s-Cs,{b-d-c-^s)Pt-sf)ds- [\^s,c-{^s-Cs)Pt-sf)ds. 

Jo Jo 

According to Lemma |3.5| and Assumptions [Ml we can find a constant Ci such that 
{b — d — c ■ ^t)Pt-sf /Cl and c/Ci belong to x {u, u}). We deduce that for all 

t 

1(6-0, V^)l < Wl(Co, Co)+C'l f sup Wl{^r,Cr)ds+Cl sup (6,1) [ sup Wl{^rXr)ds. 

Jo rg[o,s] '’sIoTjq] Jo rS[ 0 ,d 

_ (4.20) 

Using Gronwall’s Lemma and (4.19), we conclude that there exists a constant C{T^^) 
such that 

sup Wi(6,0) <c(M)Wi( 6 ,Co). 

^e[o,T 5 j 

Choosing = 6/C{J^ we fi nd for all Co with >Vi(Co, Co) < a^o, sup^gp^r^J >Vi(Cr, Cr) < 

5. In addition with ( 4.18| ) and (4.16), we find that for all Co with W’i(Co,Co) < Q^^o, 

sup sup >Vi(Ct,t '0 > 
ie[0,T5(,] v&A{To) 

Remark that adh{{B\ x R(0, e)) is a compact set of Mp{X x {tt,u}) as 

X is bounded. It is covered by p)xB(o,e))-®(Co, ct^o)- extract a finite 

cover Finally, dehning Tq = maxj=i.,„, we conclude : for any Co with 

>Vi(C 66 ) e]/3,7[ and (Cg,!) < e, we have 


inf sup >Vi(Ct,t't) > d- 
^^XTo)te[ 0 ,To\ 


Proof of Lemma 4.3 is now complete. 


Proof of Lemma 4-4 Lemma [4.3| gives Ti such that 


limsup^log sup P^K(r A A > Ti) < —1. 

K^+co K i/^e(B(C“,7)xS(0,e))nM|f ° 


□ 


(4.21) 


Thus we limit our study to the time interval [0, Ti j. S ince the considered initial states 
satisfy {i'q , 1) < ((C“, 1) + 27 + 2e), using Lemma [ t^ we hnd A" > 0 such that 

limsup—log sup ®ap (t'/^, 1) > A) < —1. (4.22) 

K ^+00 P i/fe(S(C“,7)xB(0,e))nM|f ° t£[0,Ti] 
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Let M > ((^“, 1) + 27 + 2e) V N and 

A = {v e !]){[{), Ti],Mf{X X {n,n})) | e [OjTi], < 7, sup < e 

te[o,Ti] 

and sup (4.23) 

ie[0,Ti] 

For all Vq , and K, 


> 7,T < Tf A 5f) 


K 


yK\ 


< F,k(t < Ti, sup ,1) > M)+F,k{t < E A)+F,k{t >Ti,t< Tf A 5f) 


K 




^K\ 


te[o,ri] 


<Px( sup {u^A)>M)+F.{v^^A)+F.{T^Tt^S^>T^). 

° te[o,ri] “ “ 

Then, we use Theorem 3.1 the dehnition of A (4.23), (4.21), and ( 4.22[ ) to find 


1 


limsup — log sup F^k{v!^ E < max 


(-1,- inf 
{ ugAfo&C } 


X—>+00 A 

where C = adh(B{^'^,3p) ^ -®(0t))- There remains to prove that the r.h.s. is strictly 
negative. Proposition 4.2 implies that any solution {^t)t>o to (2.1) with ^0 £ C satisfies 
the property : 

if sup 1) < e' = 2e, then T) < 7/2- 

te[o,ri] 


We deduce immediately that for any v ^ A, with uq E C, if ^ is the solution to (2.1) with 
^0 := 

sup WiiutAt) > sup max{>Vi(z.^“,Cr),>Vi(t'r,?D} > ? A e. 


ie[o,ri] 


te[o,Ti] 


We conclude the proof using Proposition 3.4 and a compactness argument similar to the 
one at the end of the previous proof. □ 

Proof of Lemma 4-5 Let us fix / E and study the following difference, using 

the construction of the process 


i(T“./> - A-'‘,/)i < f 


if(xi) - f(x‘)i + Y2 i/W)i+ E i/« 


1 ^^^notdead 


jgjydead 


jgjYborn 


where TVnotdead jg gg|. q£ indices of the individuals with trait u alive at time 0 and not 
dead during [0,t]; represents the set of indices of the individuals with trait u alive 

at time 0 and dead during [0, t]; and is the set of indices of the individuals with 

trait u born during ]0,t]. As / E C^^P{X), 




E 




Ardead i Arborn 
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As this is true for aii / G C^^P{X), we find the same upper bound for , 1 ^^’^). 

We use now the stopping time 

'Tn = inf{t > 0, 1)1 > N}. 

On {t^ > t}, is stochasticaiiy bounded by ^ \^t ~ ^h\ + 

where {(W))4>o}ig{i,.xAr} ^-re KN independent reflected diffusion processes driven by 
Equation (1.2) with the diffusion coefficient and {V{t))t>o is a Poisson process with 
intensity (b + d + Nc)KN. Finaiiy, 


( sup > p) 

ne[o,r] ^ 


< < 


T) > T, sup > p), 

^ te[o,T] ^ 


KN 


<nr^ <T)+¥i^Yl \Xl-Xi\>P-)+¥( sup ® 

^te[o,T] He[o,T] 


K 


(4.24) 


— r» / ’ 


Using Lemma 3.2, we can fix G N such that 


1 


iimsup—iogP(rj() < T) < —C. 

K^+oo X 

Let us now consider the second term of (4.24). For x = {x^,x^^) G , we denote 
the probabiiity under which (A^q, .., ) is equai to x by P^. Let T be the stopping 

time T = inf{s > 0, \Xl — Ag| > pK/2}. Using the Markov property, we And 


KN 


pK\ 


fJ sup =Pj.(T<r) 

Ue[o,T] ^ ^ ^ 


KN 


pK\ 


KN 


<IPx(T<r, ^ 

i=0 


T <r,^|A^-A^l > 


2 = 1 


KN 




i=0 


pN\ 

4 ) 


KN 


+ lEj 


'-40 


■ W 


El^i 


2=1 


T-s-Xq\ > —j 


4 / 
Ixeds 


< 2 sup Fy 




i=0 


(4.25) 

The aim is thus to And an upper bound on the last term for any x G X^^ and any 
s G [0,T]. Using Markov’s inequality, 


KN j. ^ KN 


i=0 


i=l 


e 


(4.26) 
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If we denote the kernel of the semigroup of the reflected diffusion process by y), 
Part 3 in |32] and the fact that X is compact imply that there exist two positive constants 
Ci,C 2 such that for any x,y € X, 


Cl 

p'^{x,y) < -^e , 


(4.27) 


Thus, using (4.27) then a change of variables, we find that there exists 6*3 > 0 indepen¬ 
dent from s such that 


E., 


e 


< 


/ 

JR" 


Cl \^-vr \v-A 


.M 


„d/2 


g 02 ® e dy < / Cie *^2 = C 3 < -koo. (4.28) 


We deduce with the last line in (4.25), (4.26) and (4.28) that 
P.f sup g \Xt - XJI > 


te[o,T] 


where — ln(C 3 )A^ tends to -|-oo when T tends to 0, i.e. there exists Ti such that for 
alir<ri, ^ - ln(C 4 )iV > c. 

Finally, concerning the third term of (4.24), let r = log 


2{h+d+cN)NT 


pf sup 

He[o,T] 


m 

K 



< f,{-K{eX-(b+d+cN)NT{e--l))) ^ 

< e(-^(f h( 2(.+.4iV)iVT )-l]+(^+'^+^~^)^^)) . 


There exists T 2 < Ti such that for all T < T 2 , f 


log 


cN)NT > C. We conclude the proof : for all T <T 2 , 


2{b+d+cN)NT 


limsup — logP( sup >Vi(z 2 ^’“, t',f'’“) > p) < —C. 
K^+00 te[o,r] 


-|- (6 -|- d -|- 


□ 


5 Survival probability for a branching diffusion process 

In this part, we study a distinct model which is a branching diffusion process. This model 
is correlated with the other one as explained at the end of Section Any individual is 
characterized by its location XI £ X which is described as before by a diffusion process 
normally reflected at the boundary of A, ( 1 . 2 ), with the diffusion coefficient m > 0. 
Moreover, each individual with location x £ X gives birth to a new individual at rate 
b{x) and dies at rate d(x). Those rates are assumed to have some regularity that we 
detail. 
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Assumption 5.1. 6, d are two Lipschitz functions, b is a positive function and there 
exists b, b, d such that for all x ^ X, h < b{x) < b and 0 < d{x) < d. Moreover, d is a 
non-zero function. 

Let Mt denote the number of individuals at time t. We describe the dynamics of the 
diffusion process at each time by the finite measure 

Mt 

i=l 


The aim of this part is to describe the survival probability of the population using its 
parameters. Let 

To = inf{t >0,Mt = 0}. 

The first Theorem concerns the survival probability of the population assuming that, 
initially, there is only one individual at location x. This probability is characterized as 
a solution to an elliptic differential equation on the location space X. Remark that the 
location of the first individual plays a main role. Indeed, if the first individual appears 
in a place where the growth is low or negative, it has a high probability to die with no 
descendants. 

We denote the probability measure under which r]o = dx by 


Theorem 5.2. Let H be the principal eigenvalue of the elliptic operator mAx- + (6 — d). 
with Neumann boundary conditions on X, see (2.2). If H >0, there exists a unique 
positive C'^ -solution cf* to the elliptic equation 


0 = mAx4>*{x) + (6(x) — d{x))4>*{x) — b{x)4>*(x)'^ , Vx G X, 
dn4>*{x) = 0, Vx G dX, 


and 4>*{x) = limt_^oo IP5^; ["I'd ^ t] for all x £ X. 


If H < 0, (5.1) has no non negative solution, we set (jf = 0, and limt_>.oo [Tq > t] = 0 
for all X G A. 

The second result of this part estimates the probability that the population size is of 
order K after a logarithmic time logiL. For all e > 0 and K G N*, we set 


Tex = inf{t > 0, (? 7 f, 1) > eK]. 

Theorem 5.3. Lete > 0 and {tK)K>o be a sequence of times such that liuiK^+oo / ^og{K) 
+00. Then for all x € X, 


lim P^JTei^ < tx] = 
K^+00 

The end of this part is devoted to the two proofs. 
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Proof of Theorem 5.2 We first study the behavior of the probability P 5 ^[To < t]. We 
denote the time of the first event (birth or death) of the population by Ei. The law 


of El is given by [^i < t] = 


foidiX,) + 


. We set 


I(s) := f^(b(Xr) + d{Xr))dr. Using the Markov property of rj, we obtain 

E5,[To<t] = Es^ tEi<t^{ME^=0} +'^Ei<tT^{ME^=2}^2Sxj^^[^Mt-Ej^=o] 

(d{Xs) + b{X,)Es^^ [To < t - s]2) 


= E, 


L Jo 


where X under P^; is solution to (1.2) with the initial condition x and the diffusion 
coefficient m. Thus g{x,t) = P 5 ^[To < t] satisfies for all x £ X, and all t > 0, 


g{x,t) = Ej 


Ljo 


..) + b{Xs)g{Xs,t-s)^)e-^^^Us 


, V(x, t) £ X X M"*", 


(5.2) 


g{x, 0) = 0, Va; G X. 


Using Gronwall’s Lemma for bounded functions, we deduce immediately that (5.2) has 
a unique bounded solution. 

We will now show that there exists a unique C^-solution to 

dtf{x,t) = 'm^xf{x,t) ~ {b{x) + d{x)) f{x,t) + d{x) + b{x)f{x,t)‘^, V(x,t) £ X x M"'', 
< dnf{x, t) = 0, V(x, t) £ dX X M’*' 

_ /(x,0) = 0, Vx G T, 

(5.3) 

such that / G x M"*"), is positive and is smaller than 1 by using super- and sub¬ 

solutions arguments. Indeed, let F(x, /) = — (6(x) -|- d{x))f + d{x) + b{x)f‘^. We easily 
see that / = 0 and / = 1 satisfy: 

f dtl < mAxl + E{x, 1), V(x,t) G T x IR+, 

\ dtf > mA^f + E{x, f), V(x,t) G T x M+, 
l{x,0) < /(x,0) < /(x,0),x G T, 
dnf{x, t) <0 < dnf{x, 0), X G dX, t £ M+. 


That is, / (resp. /) is a sub-solution (resp. super-solution) to (5.3). Moreover, E and 
dfE belong to C'(T’ x M) and T is a Lipschitz function with respect to x by means of 
Assumptions 5.1 We apply Theorem 4 of Chapter III of |29) to deduce that (5.3) admits 
a solution / G (7^3(^ x M+) satisfying 0 < / < 1. The uniqueness of the solution is a 
consequence of the maximum principle. 

The next step is to use a Feynman-Kac formula to deduce that / is also a solution 


to (5.2). Indeed, for X solution to ( |1.2| ), let, for all t > 0, and s £ [0, t], 

H{s,Xs) = f{XsX-s)e-^^^\ 
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Applying Ito Formula to H{s,Xs), and using (1.2), (5.3) and the fact that dnf{x,t) = 0 
for all X G dX, we find for all s G [0, t], 


H{s, Xs) = H{0, Xo) - (d(X,) + b{X^)f{X^, t - af) 


+ / V2m{dxf{X^,t-(^))e 
Jo 


(5.4) 


We take the expectation of (5.4) for s < t. The expectation of the last term is equal 
to 0 as (f^ dxf{X„,t — is a martingale. In addition, Kx[H{0,Xq)] = 

f{x,t)- It stays to make s tend to t using the dominate convergence Theorem. As 
Ea;[Fr(s, Xs)] —)> Xt)] = 0, we deduce that / is a solution to (|5.2|). (5.2) admits 


s — 


a unique bounded solution, thus both solutions are equal, i.e. < t] = fix,t)- 

Finally, we deduce the survival probability limt_^+oo(l — f{x,t)) using results on 


Equation (5.3) described in |2] and in Theorems 9 and 11 of Chapter III in |29]. Indeed, 
they prove that if Ff >0, there exists a unique positive solution (p* to the elliptic equation 


(5.1), and that 4>{t,x) = 1 — f{x,t) = > t) tends to (f)*{x) in C^{X) as t —>• Too. 


Moreover, if FF <0, the unique solution to (5.1) is the zero function and 4’{t,x) —)• 0 
uniformly in A as f —>■ +oo. □ 

Proof of Theorem \5.3[ First, we split the studied probability into three parts : 

< tx) = < tA',loglog(FF) < To < Too) ^ 

+ < tx, loglog(FF) > To) + P 5 „,(Teii- < tx, To = Tec). 

Let us start with the first term of (|5.5|) : 


start with the first term of ( |5.5[ ) : 

{T^x < tx,loglog{K) < To < Too) < P 5 ^(loglog(FF) < To < Too) 


—r 

K^-\-oo 


0 . 


The second term of (5.5) will be treated using a comparison with a pure birth process. 


Let us consider a birth process with constant birth rate b and started with only one 
individual. T^x denote the first time when N, the population size of the process, is 
greater than eK. 

n.(T.x < ti^,loglog(FF) > To) < FsA'^eK < loglog(FL)) 

<Pi(t,i^ <loglog(FF)) 

< IPi(-^iogiog_ft: > eK) 

< ^ 0 . 

K^-\-oo 


There remains to deal with the third term in (5.5). Note that if Ff < 0, Theorem (5.2) 


implies that the third term is equal to zero, and the proof is done. 

From this point forward, we assume that FF > 0. Let F be a positive eigenvector of the 
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operator mAx- + {b — d). with Neumann boundary conditions on dX associated with the 
eigenvalue H. Thanks to Ito’s Formula, we find 


pt pt 

{rjt,e~^^h) = {r]o,h)+ / {r]s,e~^''mAxh+{b-d)he~^''-Hhe~^'')ds+ / xh{Xl)e~^‘'dBl, 

Jo Jo 


As mAxh + (b — d)h = Hh and Xxh is bounded on X, h))t>o is a martingale. 

Moreover it is positive, so, it converges a.s. to a non-negative random variable that will 
be denoted by W. Obviously, {Tq < -|-oo} C {W = 0 }. Our aim is to prove that this is 
an a.s. equality. As done in the previous proof, we denote the time of the first event of the 
population by Ei and we set I{s) := jQ{b{Xr) -|- d{Xr))dr. Using the Markov property 
and the independence between individuals, we find an equation satisfied by = 0 ]: 


.\W = 0 ]=Ex 


r*+oo 


L^o 


(d{Xs) + b{X,)¥s^^ [VF = 0]2)ds 


( 5 . 6 ) 


g{x) = E,, 


■ r+oo 

/ e-^^%{Xs)giX,)i 2 -giX,))ds 

Jo 


,Vx G A. 


( 5 . 7 ) 


Finally, g{x) := E5„,[VF > 0 ] = 1 — = 0 ] is solution to 

f+OO 

UO 

Let us show that there exists at most one non-zero solution with value in [ 0 , 1 ] to Equation 
( 5 . 7 ). Let gi and g2 be two such solutions. We define 

7 = sup{7 > 0,51 (x) - 752(3;) > 0 ,Vx e X}. 


Assume first that 7 < 1. 


5i(x)-752(x) = Ej 


r-\-oo 


e-"W6(A,)[2(5i(X,)-752(X,))-(5i(X,)2-752(X,)2) 


L30 


ds 

( 5 . 8 ) 


As 7 < 1, 5? - 75^ < (51 - 7fi'2)(5i + 752) and we find 

r*+oo 


51(3;) - 752(3;) > E,j 


r+oo 

/ e-^(^)6(A,)[5i(X,) -752(X,)] [2 - 5 i(X,) - 752(X,)] ds 

Jo 


Moreover, by the definition of 7, there exists xo £ X such that 51 (xq) — 752 ( 3 ;o) = 0 , so 
/•+00 


0 = E 


Xo 


I JO 

Thus for a.e. s G IE 


e-A-)6(X,) [51 (A,) - 752(W,)] [2 - 51 (A,) - 752 (A,)]ds 


, Exo-a.s^, b{Xs)[gi{Xs) - 752(As)][2 - 51 (A^) - 752(^5)] = 0. 
Let us note that for all x G A, for i = 1 , 2 , (| 5 . 6 |) implies that 


1 - giix) > Ej 


1 


,-I{s 


U{Xs)ds 


= [Mei = 0] > 0, 
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that is, for all x G Af, 2 — gi(x) — jg 2 (x} >0. As 6 is positive, for a.e. s G M"*", Pa;o-a-s- 
gi(Xs) — jg 2 (X,s) = 0. In addition with the fact that jg^ — gf = Agii^ ~ l) ~ (5i ~ 
192 ){gi + 752 ), implies that, 

0 = 51(3^0)-752(2^0) 

r*H-O0 


= E, 


XQ 


r-\-oo 

/ e~^'^'''^h{Xs)-ig2{Xs)[l--i\ds 

.Jo 


Using the same argument as before and that 1 — 7 > 0, we deduce that for a.e. s G M"*", 
Pa;Q-a.s., g 2 {Xs) = 0. Moreover, under ^ 3 ,^,, the random variable Xg has a density with 
respect to Lebesgue measure, that is, for Lebesgue-a.a. x G X, g 2 {x) = 0. This is a 
contradiction with the fact that g 2 is a non-zero solution. 

Finally if 7 > 1, we define instead 7 ' = sup {7 > 0,g2{x) — 751 (x) > 0,Vx G A} < 1 
and we use symmetric arguments to reach a contradiction. Thus, there is at most one 
solution to (5.7) with values in [0,1]. 


The next step is to show that (j)*, which is solution to (5.1 1 , is also solution to (5.7). Let 
us write = l — which satisfies 


0 = mA^f* - (6 + d)f* + d + b{f*f, on A, 
dnf* = 0, on dX. 


(5.9) 


We apply Ito’s Formula to f*{Xt)e fo(^i^A+d{Xr))dr^ Then taking the expectation and 


using Equation (5.9), we deduce 

rt 


r{x) =e 3 


e-^^^^[diXg) + b{Xs)r{Xsf]ds 


'0 


+ E3 




Our aim is now to let t tend to infinity. Note that I{t) > ht for all t G M"*" and that f* 
is bounded by 1. Hence, we use the dominated convergence Theorem to find those two 
convergences: 


r(W)e-^W 


< E3 


.-Ht) 


t —^-|-oo 


0 


■ /*+oo 1 r ^+00 

/ e-^^^^[d + b{f*)\Xg)ds + d)ds 

Jt J Ut 


t —^“|~CXD 


Thus, we make t tend to infinity and we find, for all x G X, 

r+oo 


rix) = E, 


/• + 00 

/ e-^^^^[d{Xg) + b{Xs)r{Xsf]ds 

Jo 


Since 4>* = 1 — f*, we conclude that (p* is a solution to (5.7). There exists at most one 
non-zero solution to (5.7), thus we have either > 0] = 4>*{x) for all x G X, or 

> 0] = 0 for all x G A. Using Ito’s Formula, it is easy to check that in the case 
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H > 0, ((??t,e ^^h))t>o is bounded in L^. So this martingale is uniformly bounded and 
it converges in to W, hence = h{x) > 0. Finally, 

> 0] = (P*{x) = PiJTo = + 00 ] ^ {PF > 0} = {To = + 00 } a.s. (5.10) 


Thus, on {To = + 00 }, 


log(eiF) ^ log((r?T,K>^e 


H > 0 a.s., 


Tf:K T^k K^ + oo 

as {r]j-^ji,h) < eiC||/i||oo, '^eK —t +00 when K tends to inhnity and PF > 0 . Hence, 

tx 


lim ;—< +00 and lim 


K ^+00 \og{eK) 


and so, the third term in (5.5) satisfies 


IF’< 53 ,C^eif < tx, To = + 00 ) = 


x ^+00 log{eK) 


tx 


= +C)0, 


log(eiF) < i^^oo • 

(5.11) 


Finally, we have shown that the two hrst terms in (5.5) tend to 0, so using additionally 

( 5 ll| ) 

lim Fs,^{T^k < tx) = (To = + 00 ) = (j)*{x). 
x^+00 


That ends the proof for H > 0. 


□ 


6 


Proof of Theorem 


2.5 


This section is devoted to prove Theorem 2.5 The structure of the proof of Theorem 2.5 


is similar to the one of [3], thus, we do not repeat all the details but only the points that 
are different. 

The hrst proposition concerns the behavior of the hrst time of a mutation , when the 
initial state is a monomorphic population. 


Proposition 6.1. Suppose that Assumptions m and (|2.4[ ) hold. Let u £ U and C“ a 
eompact subset of Mp{X x {u}) such that 0 0 C“. If Vq G C“ U Mp{X), then 


for any 7 > 0 , lim > logK, sup ,f,^6u) >7 = 0 . 

° V te [log K,sf] J 


• Furthermore, ^1™ > T^') ~ (y~ ^ fx b'^ix)p^{x)f,^{dx)^. 

Proposition |6.1| can be proved using similar arguments as those of the proof of Lemma 
2 in [ 3 . It is a consequence of the following lemma. 
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Lemma 6.2. For any a > 0, there exists > 0 such that for any ^ for any 
t > To,, < a, where {f,t) t>o is the solution to Equation ( |2.1[ ) with initial state 

f,o- 


Proof On the one hand, Theorem 1.4 in |22) implies that the density of is a stable 
monomorphic eqnilibrinm for the L^-distance. Using a proof in three steps as that of 


Proposition 4.2 we prove a Wi-stability : there exists a' snch that for any a') 


and for any t > 0 , < a. 

On the other hand, for any converges towards There exists Tio snch that 

< of/2. Using Lemma 3.5 and argnments similar to (4.5) and (4.6), we 


show that for any f > 0 , any 


sup >Vi(Ct, 6 ) < C'(i)>Vi(Co,^o), 

re[0,i] 


where {Ct)t>o is the solution to (2.1) with initial state Co £ Mp{X x {n}). Consequently, 
there exists > 0 such that for any Co £ < ol Thus, for 

any Co such that >Vi(Co,Co) < ^ny t > >Vi(Ct,C“) < «■ 

Finally, as is a compact set, there exists a hnite number of balls such that C“ C 
). Dehning Tq, = maxj=i,,„ T^i, we deduce the lemma. □ 


Proof of Proposition 6.1 First, remark that the hrst probability of Proposition 6.1 


IS 


non-increasing with 7 . Thus, it is sufficient to prove the property for any small 7 > 0. 


Let us assume that 7 satishes the assumptions of Theorem 4.1 Theorem 4.1 in the 
monomorphic case implies that : there exist 7 ' > 0 , U > 0 , such that 


sup P [R^ <SfA e^^) -A 0 , 


is dehned by (4.1). We set 2a = 7 ', then Lemma 


6.2 


sup > 2 a) ^ 


and Theorem 

0 . 


3.1 


( 6 . 1 ) 

imply that 

( 6 . 2 ) 


Using the Markov property, we deduce if K is sufficiently large such that log(iF) > Tq, 


X > log K, sup Wi (i/f, C“) > 7 


te[iogi^,sf] 


<P,x(Wi(z.£,r)>7')+E.T 


< ,C“) > i)+^uK < 5'f) 

,K Pu\ 


>iog(ii:),Wi(^'|^ ,C“)<7'} sup Wi(r'i ,C“)>7 

t&[T^,S^+T^] 

,K ^ cK 


< P X(Wi(i.^,C“) > 2 a) + sup P,o {R^ <Sf A e^^)+ P.^(e ^^ < R^ < S^) 

uo&B{tn'r 
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The first term and the second term tend to 0 when K tends to +oo according to (6.2) 
and (6.1) respectively. There remains the third term. On {t < }, the number of 

mutations is stochastically bounded from below by a Poisson process with parameter 
KqK{{hy,n-l\\h>\\L^p) which is positive if 7 is small enough. We conclude with the 
fact that < F{M.f,Kv = 0) —)■ 0, under Assumption (|2.4|). That 

ends the proof of the hrst point. 


The second point of Proposition |6.1| is easily deduced from this hrst point and Lemma 2 
in [3]. □ 

The second proposition studies the process with a dimorphic initial state. Let us 
dehne : 


6q is the hrst time when the population becomes monomorphic again, 
Vo is the phenotypic trait of the population at time 9q. 


Then, 


Proposition 6.3. Suppose Assumptions\l . 1\\2.^ and (2.4) hold, and that the initial state 


Vq G Mf{X X {u, t>}) is such that converges weakly to in Mf{X) and 


lim P(Vb = u) = 1 — lim P(Vb = u) = i?!)''“(xo). 

K^oo K^oo 


Moreover, 


Vr? > 0, lim I 

K^oo 


^o<S, 


K 


A and\f-f > 0, < 7 ) = 1. 


Proof. We set 7 > 0 small enough to use Theorem 4.1 : there exist 7 ', e, P such that 


sup 




.0 

where R^, have been dehned by (4.1), (4.3). 


{Ri^ > A Sf A -A 0, 

’^0 ~ ^ ^ R'-S.OO 


Let assume that e < 7 and that K is large enough such that qx < 7 . Then on {t < 
R^^ A A }, the process ’)t>o is stochastically bounded: 


K 


K 


Z^iT-f ^ Z^^P are two branching diffusion processes starting with one individual at location 
xq , their birth rates are respectively b'^{x){l — 7 ) and b'’{x), and their death rates are 
dA{x) + + 2c7 and d^{x) + — C 7 . 

Let us set 


= inf{t > 0, 1) > Ke}, 

Tq^p = inf{t > 0 , {Z^^P, 1 ) = 0 }, 
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and respectively associated with . 

Using same kind of computations as in Lemma 3 in [3], we deduce that 


K [9o<S^ A 




> 


I r^SUp 


< 


V 


— 0 — 


qxK 


A ) - 


T] 


qkK 


> sf 1 -: 


qRK 


A Sf A rf > 


(6.3) 

Theorem 4.1 implies that the last term tends to 0 under Assumption (2.4). The second 
term tends also to 0 as the number of individuals is stochastically bounded from above by 
a birth and death process with birth rate b and competition rate c, thus Lemma 2 in 
implies that for any <5 > 0 , there exists rj > 0 such that limsup^_^_|_oo P(5f < 

Thus the main difficulty is to evaluate the first term. Theorem 5.2 implies that 


^“0 * Ac <+“) = !- '#"'”"(*0), (6.4) 

where is the solution to the following elliptic equation on X with Neumann bound¬ 
ary condition 

+ [b^x) - (Fix) - • e + - b^{x)^'^’^^{xf = 0 . 

Let us show that this solution is close to Theorem |5.2| implies that is positive 
if and only if 

ffU^VU 

H--c^-.e + n = H^-^^ + ci>0. (6.5) 


First case: < 0 (Point 1 in Assumption 2.3) 

We can find 7 small enough such that (6.5) is not satisfied, thus = 0. In 

addition with (6.3) and (6.4), we deduce that 

V 


K{eo< A ^,Uo = u,W,{ug,C) < 7 


A'^+oo 


Proposition 6.3 is proved for this case. 


Second case: > 0 (Point 2 in Assumption 2.3) 


Hence (6.5) is satisfied for all 7 > 0. Let C be 
{b‘^{x) — (F{x) — -|- C7)/ — b^{x)f 


yex b" 

^ We have 


C7 




and set C^{f) = nf Axf + 


> 0 , 

£^((1 + (£)<(.’'“) = (1 + C )<^™[ c 7 - C 6 >™] < 0 . 

As is the unique solution to C'^{f) = 0, we deduce the following inequalities from 
a comparison theorem (see for example Theorem III.5 in |29| 1 : for any x € X, 


{1 + 0 > <^'^’^“(x) > (^"“(x). 


( 6 . 6 ) 
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We split the end of the proof into three steps regarding as the proof of Lemma 3 in |3] . 


Let ns fix a > 0. (6.3), (6.4) and ( 6 . 6 ) imply that if K is large enongh, 

7] 


..K I < 


0 - ^ A A i?:;' A ) > 1 - </."“(xo) - «. 


On the other hand, if K is large enongh : 
V 


V, 




A Sf A 


> 


^2^0 \ Ke — 


1 

qxK 




dnf 


> 1 - 


qkK 


A 5 f A > R, 


QkK 


K 


e — 


> - a. 


Tims, the n-process reaches a non-negligible size Ke, before A A R^, with 
a probability that tends to i^^^(xo)- 

Once the mntant popnlation has reached a non-negligible size, we can compare the 
stochastic process and the deterministic limiting process. Under point 2 in Assnmp- 
there exist T > 0 and 72 > 0 snch that for any .^0 £ CLdh{B, •j) x (i?( 0 , 2 e) \ 


2.3 


tion 


B{0, e)), for any t >T, 


6Gi?(0,72)xi?(f,72), 


where ^ is the solntion to (2.1 ) with a dimorphic initial state ^0 £ Mf{X x {n, v}). This 
can be proved nsing similar argnments than those of the proof of Lemma |6.2[ 


Moreover, nsing Theorem 3.1 and Proposition 3.4 


supP K sup < 72 

° \t&[o,T] ° , 


iC^+oo 


with the solution to (2.1) with initial state Vq . The two previous results and the 
Markov property imply that, if K is large enough, 

I ^ A Sf A <, 5 f > ^ + T, E . 6 ( 0 ,272), E , 272) 


■'o V ' “ qxK 
> (?!)™(xo) — 2a. 


QkK 


Finally, we use the Markov property at time + T and we conclude as in Lemma 3 
in |3]. If 72 is sufficiently small, we prove that, with a probability that tends to 1 , after 
time + T, the n-population process will become extinct before its size reaches 
the threshold and before the n-process moves away from a neighborhood of the 
equilibrium . 

That concludes the proof of Proposition 6.3 for the second case. □ 


Theorem 2.5 is deduced from Propositions 6.1 and 6.3 in a similar way to Theorem 
1 in [3] by using the transition probabilities of the jump process {At)t>o- 
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7 Numerics 


To end this paper, let us illustrate Theorem 2.5 with a numerical example. We consider 
here a set of parameters similar to the one in m and |22) . The location space is X = (0,1) 
and the trait space is lA = [0,1]. For all trait n, we consider that the growth rate is 
maximal for the location x = u. For instance, the location space state can represent 
a variation of resources, as seed size for some birds, and so two populations with two 
different traits are not best-adapted to same resources. For some bird species, a gradual 
variation of seed size can determine a gradual variation of beak size [19]. Moreover, the 
maximum value of the growth rate function is the same for all traits but when u increases, 
the birth rate function goes faster to 0, as follows : 

b{x, u) = max{4 — 160 • u{x — 0}, d{x, u) = 1. 


That is, the birds with a small trait value are more generalists than the ones with a 
large trait value : they are adapted to a larger set of seed size m All individuals move 
with the same diffusion coefficient = 0.003. The competition kernel is a constant 
c = 10 and the mutation kernel k{x, u, w) is the probability density of a Gaussian random 
variable with mean u and standard deviation 0.05 conditioned on staying in U. The initial 
population is composed of K individuals at location x = 0.5 and with trait u = 0.6. 
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Figure 1: Simulations with K = 100000, qx = W 

In Figure]^ we observe the evolution of the trait associated with a change of spatial 
niches and spatial patterns over time. After a short time (a), the spatial distribution of 
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the monomorphic population with trait u = 0.6 stabilizes. Then, in Figures (b), (c) and 
(d), we observe a phenomenon of invasion and replacement : some individuals with trait 
u = 0.515 appear, invade and finally replace the previous population with trait u = 0.6. 
Remark the change of spatial niche, see Figure (c). The locations of the population with 
trait u = 0.515 are slightly smaller than the one of trait u = 0.6. Three other phenomena 
of invasion and replacement with a displacement of the spatial niche are detected until 
the time t = 2330 (Figure (e)). In a second phase, the population evolve to become more 
and more generalists (Figure (f)) : the length of the spatial niche is increasing at each 
event of invasion and replacement. 

The simulations are computed using the algorithm presented in Parts 3 and 6 in Cham- 
pagnat, Meleard |5]. It is an iterative construction, which gives an effective algorithm 
of the process. The diffusion motion is simulated using an Euler scheme for reflected 
diffusion process. 
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